Does long-range antiferromagnetism help or inhibit superconductivity? 
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We analyze the possible existence of a superconducting state in a background with long-range 
antiferromagnetism. We consider a generalized Hubbard model with nearest-neighbor correlated 
hopping in a square lattice. Near half filling, the model exhibits a d-wave-Bardeen-Cooper-Schrieffer 
(BCS) solution in the paramagnetic state. The superconducting solution would be enhanced by 
the antiferromagnetic background if the contribution of triplet pairs with d-wave symmetry and 
total momentum (tt, tt) could be neglected. However, we find that due to their contribution, the 
coexistence of superconductivity and long-range antiferromagnetism is ruled out for large values of 
the Coulomb repulsion U. Spin-density wave fluctuations (SDWF) do not change this result. 
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1. Introduction. 

Eleven years after the discovery of the high-Tc 
superconductivityu several questions concerning the 
pairing mechanism and other important features of the 
phase-diagram of the superconducting cuprates remain 
without a satisfactory answer. Fortunately, there are 
also some aspects of the complex nature of these mate- 
rials that became quite clear. It is widely accepted that 
the antiferromagnetic correlations play an important role 
in the physics of the superconducting phase and there is 
a good amount of evidences in favor of the d-wave sym- 
metry of the superconducting order parameter. 

The proximity between the superconducting and an- 
tiferromagnetic phases in the phase diagram of the 
cuprates inspired the theoreticians of this field in differ- 
ent ways. In some theories an antiferromagnetic back- 
ground with long-range, osder is assumed to providerthe 
scenario for the pairingpO. In other class of theoriesQu, 
the pairing takes place in a background of short-range 
antiferromagnetic correlations with coherence lengths of 
only a few lattice sites. Furthermore, the different phases 
of the cuprates have been considered to be the result 
of the-competition between these two different kinds of 
ordertj. The coexistence between antiferromagnetism 
and superconductivity has been discussed in otherj^jS^ 
tems with conventional electron-|phonon interactiorO'EHI 
or in heavy fermion compoundstS. In particular, some 
Chevrel phase compounds containing rare earths, exhibit 
an anomaly in the dependence of the upper critical field 
with temperature- as the system becomes antiferromag- 
neticaly orderedtil. Thia. behavior has been explained 
using Eliashberg theoryllj and is due to the change in 
the quasiparticle-phonon interaction in the magnetically 
ordered phase. 

The model we consider is defined by the Hamiltonian 

i i 

~ X! (^o-Cjs- + h.c){tAA (1 - nic){l - rija) 

<ij>(7 

+ tAB - «ja) + rij-ff (1 - riicr)]}, (1) 

where < ij > denotes nearest-neighbor positions 
of the lattice. It was obtained from a reduction 
of the three-band extended Hubbard model to de- 
scribe th&iJow-energy physics of the superconducting 
cupratestaEJ. The values of three different hopping inte- 
grals, tAA,tAB,tBB, depend on the values of the param- 
eters of the three-band Hamiltonian, i.e. on the Cu-0 
hopping tpd, the charge-transfer energy A and the on- 
site Cu Coulomb repulsion Ud- In this paper, we con- 
sider the case tAB > tAA — tBB, which corresponds to 



the limit tpd <^ A and tpd <S J7d — A. This minimal 
one-band model was useful to study several featU|t:eS|-af 
the normal state of the superconducting cupratesEJilH. 
A more accurate description jfaquires the addition of a 
next-nearest neighbor hoppin^o, which we neglect here 
for simplicity. 

Mean-field Hartree-Fock (HF) and Bardeen-Cooper- 
Schrieffer (BCS) approximations are the simplest ap- 
proaches to the study of correlated systems. Although 
these techniques are not expected, to be valid in the 
strong coupling regime {U — > oo)E3 they can give valu- 
able insight within the weak to intermediate coupling 
ones. In particular, the transition to the insulating spin- 
density-wave (SDW) phase that takes place at finite U 
for tAB < tAA — tBB is described with the SDW-HE-ap.- 
proximation with acceptable quantitative precisiont2lE3. 
We previously studied the phase diagram of (0) us- 
ing BCS and SDW-HF for thcr^ame relation of hop- 
ping parameters considered herecil. We found that the 
correlated-hopping term gives rise to an effective pair- 
ing interaction with components in the s- and d-wave 
channels. While the s-wave BCS solution is stable for 
low densities, the BCS solution with d-wave symmetry 
exists in the paramagnetic phase near half filling, within 
the same range of densities as the SDW-HF solution. 
The latter solutions are nearly degenerate at C/ = 0. 
However, for finite U the SDW solution is the one with 
the lowest energy. 

The aim of this work is to investigate the possible 
coexistence of both kinds of order near half filling. The 
density of states in the split band structure of the antifer- 
romagnetic state has a van Hove singularity at negative 
energies and it is further enhanced by the opening of an 
antiferromagnetic gap. This situation certainly enhances 
the magnitude of the BCS-gap within a range of densi- 
ties close to that corresponding to the optimal doping 
of the cuprates. Such a favorable situation is the main 
ingredient of the first class of theories above mentioned. 
However, due to the broken symmetry of the antiferro- 
magnetic state, not only singlet but also triplet pairs are 
coupled by the attractive interaction. We show that for 
certain kind of interactions, like the one considered here, 
or the exchange interaction t/X]<ij> ' of t — J-like 
models, both contributions have opposite sign and tend 
to cancel in the strong coupling limit. We also show 
that this tendency is not modified by spin-density wave 
fluctuations (SDWF) treated within the random-phase 
approximation (RPA). This kind of behavior might be 
related to the physics of the underdoped materials. 

In section 2 we describe the mean-field picture and 
discuss the results. Section 3 contains the analysis of 
the RPA fluctuations. We summarize and interpret our 
results, and compare them with other theories in section 
4. 



2 



2. Mean-field solution. 

It is convenient to separate the correlated hopping terms 
of the Hamiltonian (nl) in one- two- and three-body con- 
tributions. The complete Hamiltonian reads 



i 

+ h {riia + nj„) + tg Tii^Tij^}, (2) 



E 

<ij>C7 



where t = Iaa, ^2 = tAA^tAB and = 2tAB-tAA-tBB- 
Details about the mean- field -decoupling of the three- 
body terms can be found iiiE3"E3. For simplicity, in what 
follows, we do not take into account_BCS terms with s- 
wave symmetry. As shown in RefEJ, these terms are 
relevant only for very high doping. 

The effective one-body Hamiltonian can be written as: 



HmF = C - ^[fief 



{Um/2 + 4t3T)]n^a 



-teff J2 ('^Lcja + h.c.) 

<ij>cr 



i,(5=x,y) 



4+5i4^ + h.c), (3) 



where 



^ef = fJ-- {Un/2 + 8t2T + UsTu) 
teff = t- t2n + <3[3r2 _^ ^2^ _ ^^2 _ 

b,s = -2t3Tipis 
C/L = -C/(n2 - m'^)/4 - 8t2nT 



4<3[4t(^2 



(4) 



The staggered magnetization (ni|) — (n^) = me'*^'^', 

= and (p^s = {cl+si:4i)^ with (5 = x,y, 

must be determined self-consistently. n = n| -I- rtj^ is the 
particle density, L is the number of sites of the lattice, 
Q = (tt, tt) and Ri denotes the atomic position. Note 
that ifix — —fiy for the solution with d-wave symmetry. 

In the paramagnetic phase to = and the anmlitude 
of the d-wave BCS gap ipix is independent of Vc3. The 
critical temperature Tc, as a function of doping is shown 
in Fig. 1. As within this range of densities the SDW 
solution (m ^ and ipi^ — 0) has lower free energy that 
the BCS one, we consider a mean field solution with both 
order parameters different from 0. Due to the symmetry 
of the SDW background, the values of fix depend on 
whether i belongs to the sublattice in which the T or | 
spins are the majority ones. The Fourier transforms of 
the self consistent parameters of this solution are: 



L 



k 



feT''-(fc+Q)i))] 



L 



k 



(5) 



where denotes a summation over k within the re- 
duced Brillouin zone, ip'^ is, in general, a complex quan- 
tity. There is also an imaginary term in r, which involves 
mean values of the form {c].^Q^Cka)- However, all the 
parameters are real for the lowest-energy solution. Note 
that iy9*3 is the mean value ofpUae 3^ = projection of 
a nearest-neighbor triplet pairll3. This is a consequence 
of the breaking of time-reversal symmetry in the SDW 
state. 

As usual, the normal terms of the mean-field Hamilto- 
nian are diagonalized by the canonical transformation: 

Cka = Ukl'iJ -^«fe7i^\ 



Uk 



Vk 



-(1--), 



(6) 



where 7|,g.\ T^^'' are defined on the valence and con- 
duction bands with dispersion relations — ±_Efc, re- 
spectively, where Ek = -\/e| + A^, — ~2tef{coskx + 
cosky), and A — m{U/2 + ^t^r). The anomalous terms 
are not exactly diagonalized by M). For example 



+ UkVkhi-^''l^V!-llV'^-H) 



( + )t^,(-)tN 



(7) 



Below half filling and for sufficiently large magnitude of 
the charge gap A, the interband pairs can be neglected. 
In the two-band basis (||) the mean-field Hamiltonian (^) 
results: 



ka a 

(8) 



3 



with a = +, — , and 

A^'' = (a + 2ukVk^^) 4t3r(cos - cos ky). (9) 

After a Bogoliubov transformation, the set of self- 
consistent equations can be cast as: 

" = 1-ZE ||(l-2/(A^)), 

ka ^ 

1 ^ CfcCOsfc^ Of(\'^\\ 

rr^ = m (| + 4t3r) ^ ^ i- a || (1 - 2/(A^)), 
^'-^il «(l-2/(A^))^ cosA;., 

= -2^E (l-2/(A^)) ^ cosfc,, (10) 

ka ^ 

where = ViW+Wf7 and /(A^) = 1/(1 + 
exp (/3A^ ) ) , with the temperature T = 1 //3 in units where 
the Boltzmann constant ks = I- Notice that in a SDW 
background ip'^, ip'^ depend on U through Uk,Vk and 
A^. It is apparent from Eqs. (|9|-|lO|) that these quantities 
contribute with opposite signs to the superconducting 
gap. At first sight it seems to-be vahd to neglect the 
contributions of the triplets ip^n to the BCS gap. Under 
this assumption the dependence of the mean-field criti- 
cal temperature as a function of doping is shown by the 
solid circles in Fig. 1. The maximum of Tc occurs for 
values of the chemical potential fi^f equal to the energy 
of the van Hove singularity of the density of states of 
the SDW solution. As this density of states is enhanced 
with respect to the non-interacting one, the value of the 
maximum Tc is, thus, enhanced with respect to the value 
corresponding to the paramagnetic solution for the same 
parameters. However, it can be easily seen fromEq. (^) 
that (p^ — as U ^ t. In other words, as the stag- 
gered magnetization increases to values closer to that of 
the Neel state {m — 1), the absolute value of the contri- 
butions of triplet and singlet states become equal and the 
superconducting solution disappears. In the weak cou- 
pling regime, the superconducting gap does not vanish. 
It is, however, highly reduced with respect to that of the 
paramagnetic phase. For example: for optimum doping, 
ip° = 0.226 X 10-2, pQ = -0.198 x 10"^ for U = 6t, 
while ip° = 0.122 X 10-^ V3<3 = -0.112 x 10"^ for 
U = 8t. When the contribution of triplets is neglected, 
= 0.283 x 10~i for optimum doping and U = 8t. The 
maximum value of ipx for the paramagnetic solution is 



(px = 0.187 X 10~^, independently of the magnitude of 
U. For U > 8t we were not able to find the numerical so- 
lution of p^ , pfi . One might expect that including mag- 
netic fluctuations, the singlet contribution, and with it 
the superconducting order parameter is enhanced. This 
is discussed in the next section. 



3. RPA fluctuations. 

The complete treatment of the fluctuations of the anti- 
ferromagnetic background of the doped system is a quite 
hard task. The doping on the Neel state is expected to 
introduce instabilities towards the formation jof incom- 
mensurate and spiral SDW kind of orderingEj as well 
as quantum transitions to a quantum-disordered phase 
with antiferromagnetic short-range-iorder, depending on 
the topology of the Fermi surfacecil. In this section, we 
de not consider these effects at all, and we concentrate 
in the study of the usual magnon SDWF, assuming that 
the doping does not affect the long-range order of the 
background. Even when this assumption is only strictly 
valid at half filling or for a strongly underdoped system, 
it is frequently used in several theories of the high- Tc, as 
mentioned in the introduction. For the Hubbard model 
at half filling, these SDWF within the RPA approxima- 
tion were shown to be enough to recover the value of 
the local magnetizalian {Sz) of the Heisenberg model in 
the limit of ?7 raO. These fluctuations were further 
proposed tcwBediate the pairing interaction in the SDW 
backgroundS'Ea. However, a more detailed studjji-pf the 
effective interactions for the Hubbard modelOEZI indi- 
cated that in fact they do not provide any pairing mech- 
anism. In what follows, we analyze if they can correct 
the cancelation effect between pairs with total momen- 
tum and Q presented in the previous section. In the 
following, we restrict our study to T = 0. 

The three-body terms are reduced to two-body ones in 
this level of approximation following a^^imilar procedure 
as with the Hartree-Fock decouplingE2l. It is found (up 
to a constant): 

<ij>(y 

~ {CjcrCi'jC'j^^Ci^a + h.c) 

+2 T {n,„nj„ - c\^c]_^c^-„Cj^)\. (11) 

The terms of the first line of Eq. (|ll| ) contribute only to 
the renormalization of the band- width and of the chem- 
ical potential, and for low-densities to the BCS-s-wave 
solution, but they do not couple with the SDWF. In the 
SDW background, the first term of the second line of Eq. 
( pT| ) gives the same contribution as At^r J2<ij> ^i^j 
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plus a nearest-neighbor repulsion, which is also reflected 
in the expression of the charge gap of the SDW solution. 
The second term, can be written as 



E 

<ii> 



{s+s- + srs+), 



(12) 



with SJ' — cJ|Ci^. Putting both latter contributions 
together, a Heisenberg interaction with antiferromag- 
netic coupling J = At^^T is obtained. The pairing terms 
come from Eq. ( [l^ ) and they couple with the transverse 
SDWF. Longitudinal SDWF do not contribute to the 
pairing mechanism within the present approach. Thus, 
we concentrate on the transverse channel. Writing Eq. 
( p^ ) in the two-band SDW basis (Eq. (||)) and keeping 
only the terms of the valence band, which are the rele- 
vant ones for the superconductivity bellow half filling, it 
is found 



k,q 



AUkVkUk+qVk+q) 



with 



Uq = -4:t3T {cosqx + cos qy). 



(13) 



(14) 



While the first term contains the contributions of pairs 
with 0-momentum: 



^l+Q\'^''~{k+Q)l^~{k+q+Q)l^{k+q+Q)t^ 



(15) 



the second term is obtained from the contributions of 
pairs of triplets with momentum Q: 



^h^-{k+Q)i^~(k+q+Q)lC(k+q)P 



(16) 



As discussed in the previous section, both terms tend to 
cancel each other. 

The spin susceptibility of the SDW state is a 2 x 2 ma- 
trix with diagonal elements (9i <?; '^)j (9 + Q; 9 + 

Q; ti;)EfE3~LZl. The matrix elements of the bare suscepti- 
bility with total momentum Q (off-diagonal matrix el- 
ements) do not contribute in the static case considered 
herecj. Away from half filling, the full RPA suscepti- 
bility, has interband as well as intraband contributions. 
The Goldstone modes of the Neel state are obtained from 
the poles of the interband part of the full susceptibility. 
In this case from the solution of: 



0, 



0, 



(17) 



with Uq given in Eq. ( p^ and xt~{l->^)-> the diagonal 
matrix elements of the bare susceptibility. This equa- 
tion reduces to the gap equation of the SDW mean-field 
solution, for cj = 0, q = Q, indicating the consistency 
of the approach. The intraband contribution accounts 
for the instabilities of the SDW antiferromagnetic back- 
ground with long-range orderr-sttid for the transition to 
the quantum-disordered stater iri and we do not con- 
sider them in the present work.|JVhen treated within 
the magnon-pole approximationBEZI, the interband terms 
result 



1 l-Vq 



1 



2 \l 1 + rjq u! — riq + id 

1 



X+ {q,q + Q;uj) 



1 



1 



2 uj — riq + i6 



LU + Clq — iS ' 



with Q„ 



2J, 



1 



Vq, and r]q 



(cos qx 



(18) 



'%)/2. 



The effective antiferromagnetic exchange coupling is de- 
fined from the dispersion jrelation of the spin waves. 
Following usual proceduresoc^O, it is obtained for the 
present case J^f = ^t1^/U + 4t3r. From Eq. (p^), 
it can be seen that at half filling r ~ (coskx) / ^ ~ 
2/(7r^A). Thus, in the limit U ^ 00, t^f Iab and 
Jf,f At\g/U , as expected from the results of a canon- 
ical transformation on Eq. (1) for large U . 

The fermion-fermion interactions constructed from the 
corresponding fermion-magnon ones (see Fig. 8 of Ref.S) 
are: 

1 ' 

^ E [- ^(«' + /2(fc, 9) V{q + Q,q + Q) 



k,q 



+ h{k,q) {V{q,q)-V{q + Q,q + Q))] 
J-)t^(-)t A-) (-) 



(19) 



where 



/i(fc,g) = ul^qUk 
f2ik,q) = vl+qul 

f3{k,q) = 2 Uk+qVk+qUkVk, 

V{q,q) = {U -Uqfx+-{q,q-iU 



2 2 

' '"k+q'^^k^ 

2 2 

'^k+q'^k^ 



0). 



(20) 



As in the case of the bare interaction, we kept only terms 
on the valence band. The first (second) line of Eq. ( pj| ) 
contains the contributions of pairs with total momentum 
{Q). The bare interaction ( p^ has components in the 
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d-wave- as well as in the s-wave-channel. The relevant 
terms for the c?-wave part, are those with k near (0, tt) 
and q near (0, 0) as well as symmetry-related points. 
SDWF would help to superconductivity in the case that 
the second line of Eq. (^9|) tends to cancel the second 
term of Eq. (^3|) while the first line of Eq. (1^ vanishes 
or has the same sign as the first term of Eq. (13^ for these 
wave vectors. For A » _Bfc A and fj,{k,q) — > 1/2. 
Expanding the structure factors fi{k,q), f2{k,q), the 
first line of Eq. ( p^ can be written as 

1 



-2 iViq,q)-Viq + Q,q + Q)) - 
iV{q,q) + V{q + Q,q + Q)). 



(21) 



The first term of Eq. ( |2l| ) exactly cancels the second 
line of Eq. dH). The remaining term, which comes from 
pairs with momentum 0, can be expanded for small q. 



V{q,q) + V{q + Q,q + Q) 



Jef 

~ 8t3T(4t3T + U) 

ek^k' (vk • (5k)2 



(22) 



where vl — dek/dki\i^g and 5k ^ ^k' — k — ko ^ q. 
The latter quantity is vaxushingly small except for ko ~ 
(7r/2, 7r/2), in which caseEj'Ell, the effective interaction in 
real space is a local repulsion in the triplet-channel plus 
a long-range dipolar interaction and it is expected |-tt) 
induce spiral distortions in the antiferromagnetic stateO. 

Thus, SDWF of the Neel state treated in the RPA 
approximation, do not modify the picture obtained at 
the mean-field level. 



4. Conclusions and discussion. 

Our results can be summarized as follows. We stud- 
ied an effective one-band model for the superconduct- 
ing cuprates in the underdoped regime. The model has 
an effective pairing interaction and a BCS- mean-field 
solution with d-wave symmetry in the range of doping 
< 6 < .5, which might be relevant for the pairing 
mechanism of the cuprates. However, in the mean-field 
level, this solution has always higher energy than the 
antiferromagnetic SDW one. For the case of the t — J 
model, for which the results of strona^coupling mean- 
fielcO as well as numerical techniquescj, indicate that 
it exhibits superconductivity, the instability of the Neel 
state is found also aL_very high doping within some 
mean field approachesEJ. In other more phenomenolog- 
ical approaches, antiferromagnetic long-range order has 



been used tQ|Simulate the short range antiferromagnetic 
correlationsu'EI. In this work we explored the possibility 
of coexistence of both kinds of order: antiferromagnetism 
and superconductivity. We found that neglecting the 
contribution of triplet-pairs with momentum Q as done 
in other approaches for simplicityEI, there is an enhance- 
ment of the d-wave superconducting solution due to the 
modified density of states of the antiferromagnetic state 
with long-range order. We found, however, that both 
contributions are equally important and tend to cancel 
each other. For large values of the Coulomb repulsion U 
there is no chance for superconductivity and for low to 
intermediate values the magnitude of the superconduct- 
ing gap is much weaker than that of the paramagnetic 
BCS solution. 

Our results are more transparent when analyzed in 
real space. The terms of the Hamiltonian here consid- 
ered, that cause the pairing, are the three-body ones. 
Treating them in mean-field and reducing them to two- 
body ones, a nearest- neighbor spin-flip antiferromag- 
netic interaction is found. This kind of interaction is 
expected to generate an effective attraction for nearest- 
neighbor spins in the singlet channel. In fact, in the 
strong-coupling limit, resonating valence bond (RVB) 
singlets, which are widely accepted to build up a super- 
conducting state, ace-|mainly a consequence of the same 
kind of interactionEra. The formation of singlet pairs 
is, however, very unlikely in an antiferromagnetic back- 
ground with long-range order. Spin fluctuations of the 
Neel state in the RPA level (interband fluctuations) do 
not correct this effect. For higher doping, intraband fluc- 
tuations are important and the antiferromagnetic corre- 
lations are short-range like in real materialsEJ. This sce- 
nario is better descjdbed by a nearly antiferromagnetic 
Fermi liquid pictureQ. We expect the pairing mechanism 
contained in this model to be active within such a_back- 
ground. This is confirmed by preliminary resultsEj. 

In other ibfiories, like the antiferromagnetic van Hove 
mechanisniQu or the spin-polaron modelQ, the supercon- 
ductivity is mainly determined by kinetic reasons. These 
theories assume that holes like to propagate in the long- 
range antiferromagnetic background distorting it as less 
as possible. In the first case, holes are assumed to move 
within the same sublattice and to experience an attrac- 
tive interaction ^ —0.6 J — J {Si ■ Sj) between nearest- 
neighbor holes. Note however, that the results of the pre- 
vious section and previous worksO suggest that it is not 
valid to replace the exchange interaction by a nearest- 
neighbor attraction in the Neel state. Numerical calcu- 
lations also show that the antiferromagnetic van Hove 
scenario does not represeat correctly the physics of re- 
alistic t — J-like models^. The t — J model contains, 
however, an explicit nearest-neighbor attraction of mag- 
nitude J/4, which plays a relevant role in most of the 
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strong-coupling appr In the case of the 

spin-polaron modelQ, this nearest-neighbor attraction is 
neglected and hopping of a hole between both sublat- 
tices is allowed by the occurrence of a local spin-flip in 
the slave-fermion representation used. The holes become 
paired through the hole-magnon interaction originated in 
the hopping process and the resulting Tc in this strong- 
coupling approach becomes appreciable in—contrast to 
our weak coupling results and those of Ref.EZl. It is quite 
easy to inspect the correlated hopping terms and observe 
that in the present model no terms like cLci^S'j^ are 
generated within the mean- field treatment. Such terms 
could lead to a polaron-like picture in the weak-coupling 
formalism. 

A BCS mean-field approximation does not lead to su- 
perconductivity in the ordinary Hubbard model {tab — t) 
with U > 0. There are instead, several calculations of the 
superconducting gap with other mean-field Eliashberg- 
like theories, based on the fluctuation-exchange approx- 
imation (FLEX)Ej, which suggest d-wave superconduc- 
tivity in the model. However many numerical attempts 
to find evidences of superconductivity in the pure Hub- 
bard model have been negativeE^. There are two fea- 
tures of the results based on the FLEX approximation 
that could lead one to speculate about some connection 
between these treatments for the Hubbard model and 
the results we presented here: a) within the FLEX ap- 
proximation, pairing is caused by an effective attraction 
in the c?-wave channel, which is originated by spin fluc- 
tuations in the paramagnetic phase. In the present case, 
pairing is originated by the interaction (^2|), which is 
precisely the nearest-neighbors static version of an ef- 
fective attraction mediated by the transverse spin fluc- 
tuations. Thus, within the paramagnetic phase, super- 
conductivity in the present model could have a similar 
origin as that of FLEX for thej)ure Hubbard model, b) 
As discussed in detail in Ref.E£l, the FLEX approxima- 
tion offers a quite poor treatment of antiferromagnetic 
correlations. In particular, the antiferromagnetic state 
is not recovered, even at half filling. The behavior of the 
superconducting gap within the FLEX approximatiorO 
is very similar to that shown in Fig. 1 for the BCS 
d-wave solution for Iab in the paramagnetic phase. 
Questions arise about what could occur in the pure Hub- 
bard model, in the case that the tendency towards long- 
range antiferromagnetism could be also included in some 
FLEX-like scheme. In connection with this latter point, 
it could be stressed that the effective interactions gener- 
ated by the SDWEi oLthe antiferromagnetic state of the 
Hubbard modeBcjEO, where shown not to be able to 
provide, a jjij/j effective attractive interaction in the d-wave 
channclEjE]. 
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FIGURE CAPTIONS 

Fig. 1. Mean- field critical temperature as a function 
of doping X = 1 — n for Iaa = tsB = t and Iab = 1-5^. 
Open circles correspond to the usual paramagnetic d- 
wave BCS solution. Solid circles correspond to the d- 
wave BCS solution in the SDW background for U — 8t, 
neglecting the contributions of triplets with total mo- 
mentum Q = (tt, tt). 
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FIGURE 1 



